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Outline of Paper 


¢ Overview of the validation case. 


¢ Overview of the code used for model validation: a releasable edition 
of the National Combustion Code ( also known as OpenNCC). 


¢ The results of a non-reacting case. 
¢ The results of a reacting case. 
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Workshop Case 


Computational domain and boundaries for the Volvo test case 


Bottom Wall 


0.12 m 
on Table 1. Operating conditions to be simulated: 
7 
Fuel Propane 
Oxidizer Air 
Mass Flow Rate 0.2083 kg/s 
Inlet Temperature 288.2 K 
Premixed Equivalence Ratio 0.65 
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Releasable edition of the NCC 


¢ Gaseous solver is a pre-conditioning enabled polyhedral finite—volume code. 
CFD cells can have an arbitrary number of faces and faces can have an 
arbitrary number of points. Spray solver is based on a Lagrangian scheme. 

¢ Asecond order accurate central or upwind scheme is used for spatial 
discretization of the Euler fluxes. A second order accurate central scheme is 
used for discretization of the Laplacian terms. 

¢ Temporal integration options: (1) non-iterative second order MacCormack 
scheme; (2) dual-time 2-stage modified MacCormack scheme; (3) dual-time 
3-4-5-stage Runge-Kutta scheme. 
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Releasable edition of the NCC 


¢ Subgrid models for mixing and combustion, Eulerian probability density 
function (EUPDF) and linear-eddy mixing (LEM) are available. 

¢ Four major differences between turbulence models in terms of coding are 
listed below 


Turbulence | Turbulence Eddy LG Coefficients 
Model Stresses Viscosity Destruction 
Term 


ae | Quadratic and C,,pK*/e pe RCP: Prescribed 
Cubic 
Linear C,pKo5A C.e(K)'5/A C,, C,: Prescribed 
or computed by 
LDKM scheme 
LES Linear (C,A)*p|S| N/A C,, Prescribed 


TFNS/LES Quadratic and = Min(C,pK*/e, — Max(pe, RCP: Prescribed 
Cubic (C,A)?p|S| ) C.e(K)'/A) C, C,: Prescribed 
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Non-Reacting case 


G 


f) 
YY) 
y] 


SA 


( 


4 
y 


o 


Ms 


eee 

au 
ay 

A] 


<> 
— 


a 
i 
To oa. 
wees 
‘a 
sy 


- 
#. 
fi 
UY 
Hp 


y\ 


dy 


See 


gi eines: 1) Pa 


1,667,518 polyhedrons are derived from workshop-provided 4mm-800k 
hexahedrons by a truncation technique. Its face-to-element ratio is increased 


from 3.06 to 5.43. 


TFNS 
Inlet flow rate = 0.2083 kg/s 
Back pressure = 100000 Pa 


Axial velocity, m/s 
-15.39 -0.1908 15 30.2 45.39 60.59 


Instantaneous axial velocity contours for the non-reacting case. 
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Non-Reacting case 


TFNS 
Inlet flow rate = 0.2083 kg/s 
Back pressure = 100000 Pa 


transverse (up and down) velocity, m/s 
-33.47 -19.22 4.976 9.27 23.52 37.76 
ESE  —Okh—EEE 


Instantaneous Y-velocity contours for the non-reacting case. 


TFNS 
Inlet flow rate = 0.2083 kg/s 
Back pressure = 100000 Pa 


9.839e+04 9.891e+04 9.942e+04 9.994e+04 1.005e+05 1.01e+05 


Instantaneous pressure contours for the non-reacting case. 


TFNS 
Inlet flow rate = 0.2083 kg/s 
Back pressure = 100000 Pa 


's of Z-component vorticity 
-2000 -1200 -400 400 1200 2000 


Contour: 


Instantaneous Z-vorticity contours for the non-reacting case. 
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Non-Reacting case 


Mean values of variables are ensemble averages of unsteady 
values of variables. They are computed in the code and saved 


to the disk. 


Axial velocity, m/s 
-10 2 6 14 22 30 


Mean axial velocity contours in mid-plane for the non-reacting 
case 


Pressure, PA 


9.844e+04 9.893e+04 9.94le+04 9.99e+04 1.004e+05 1.009e+05 


Mean pressure contours in mid-plane for the non-reacting case 
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#—+ Experimental 
ot TENS 


Normalized mean axial velocity 


=) 
to 


4 6 8 10 
Centerline x/D 


Centerline profiles for the normalized mean axial velocity on a 1.6m 
polyhedral grid for non-reacting case from TFNS option (i.e. u,/(16.6 m/s) ) 
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Non-Reacting case 
| 
; | 


#—e Experimental 
— _TFNS 
| 
{ i 
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| 
A 
oe 05 0 05-05 0 0505 0 0.50. 50. 
ae x/D = 0.375 x/D = 0.95 x/D = 1.53 x/D = 3.75 x/D =94 


ftaeweis *¢ fs & Ls 


D=0375 x/D=0.95 wWD=153  x/D=3.75 

Transverse profiles of normalized Transverse profiles of normalized 

mean axial velocity at different axial mean Y-velocity at different axial 
locations 


locations. 
Temporal fluctuations of variables sets for the non-reacting case are not 
post-processed due to lack of disk space to store unsteady data. 
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Reacting case 


From the previous slides, it is speculated that the reason of the poor 
agreement between the numerical results and the experimental data is 


due to the type mesh used in the non-reacting case, hence the original 
800k all hexahedral mesh from the workshop is used for the reacting 


case. 
Inlets2 2225) Outlet 
Fuel C3H8 
Air O2 and N2 
Bulk inlet velocity 17.3 m/s 
Inlet Temperature 288.2 K 
Premixed equivalence Ratio 0.65 
Exit static pressure 100000 Pa 
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Reaction Mechanism for Reacting case 


The mechanism is quite similar to the one recommended by the 
workshop. The changes are: (1) The stoichiometric coefficients of the 
reaction kinetics have been multiplied by 2 to be integers because the 
code only accepts integer stoichiometric coefficients and thus the pre- 
exponential factor of reaction constants have been divided by 2. (2) The 
global modifiers have been applied to the second reaction in order to 


keep exponential factors of concentrations the same as those 
of workshop mechanism. 


REACTIONS 
2 C3H8 


GLO / C3H8 
GLO / O02 
2CO 
GLO/CO 
GLO /O2 


cgs Cal/mol 
+702 => 6CO +8H20 1.0E+12 33000 
0.9028/ 
0.68535/ 
+102 <=> 2CO2 2.25E+10 12000 
1.0/ 
0.5/ 
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Reacting case sal 


Two turbulence models were selected for simulations, (a) time-filtered Navier- 
Stokes (denoted TFNS) approach, (b) sub-grid kinetic energy based LES approach, 
denoted K-LES. The resolution control parameter (RCP) in TFNS is set to 0.5. 
The coefficient, C,, of K-LES eddy viscosity is 0.067 initially and the coefficient, C,, 
of the destruction term of kinetic energy transport equation is 0.916 initially. The 
option to compute them (C, and C,) using the “localized dynamic kinetic energy 
model” (LDKM) was turned on later. 


Unsteady axial velocity, m/s, K-LES, 800k grid 
-15 3 21 39 57 


K-LES-LEM 


Unsteady axial velocity, M/S, K-LES-LEM, 800k grid 
21 39 57 


Instantaneous axial velocity (M/S) contours of a two-step global kinetic reaction 
model for propane 
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Reacting case 


Unsteady Y-velocity, M/S, K-LES, 800k grid 
-20 -12 4 a 


K-LES-LEM 


Unsteady Y-velocity, M/S, K-LES-LEM, 800k grid 
-20 -12 4 a 12 


Unsteady Y-velocity, M/S, TFNS, 800k grid 
-20 -12 -4 4 


Instantaneous Y-velocity (M/S) contours of a two-step global kinetic reaction 
model for propane 
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Reacting case 


Unsteady pressure, PA, K-LES, 800k grid 
9.5e+04 9.68e+04 9.86e+04 1.004e+05 1.022e+05 1.04e+05 


ES ea 
tlk es - —_ 


Unsteady pressure, PA, K-LES-LEM, 800k grid 
9.5e+04 9.68e+04 9.86e+04 1.004e+05 1.022e+05 1.04e+05 


Unsteady pressure, PA, TFNS, 800k grid 
9.5e+04 9.68e+04 9.86e+04 1.004e+05 1.022e+05 1.04e+05 


Instantaneous pressure (Pa) contours of a two-step global kinetic reaction 
model for propane 
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Reacting case 


Because a two-step global kinetic reaction model for C3H8 is used, the 
influence of the mass fraction of CO is small, the geometry of the flame 
is less wavy. 


Unsteady temperature, (K), K-LES, 800k grid 
284 627.2 970.4 1314 1657 2000 


K-LES-LEM 


Unsteady temperature, (K), K-LES-LEM, 800k grid 
284 627.2 970.4 1314 1657 2000 


Unsteady temperature, (K), TFNS, 800k grid 
284 627.2 970.4 1314 1657 2000 


Instantaneous temperature (K) contours of a two-step global kinetic reaction 
model for propane 
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Reacting case 


Experimental 


to 


Normalized axial mean velocity 


0 2 4 6 8 10 
Centerline x/D 


Centerline profiles for the normalized mean axial velocity on the 800K grid for 
reacting solutions from K-LES, K-LES-LEM and TFNS options. 
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Anisotropy(v rms/u’rms) 


Reacting case Nasal 


* Experimental 
—— K-LES 
—— K-LES-LEM 
—— TENS 


* Experimental 
—— K-LES 


Normalized Turbulence Intensity 


0 Premiet | cam ime: 
0 2 4 6 8 10 é 
Centerline x/D Centerline x/D 
Centerline profiles for the normalized Centerline profiles for the normalized 
fluctuation level of anisotropy fluctuation level of turbulence intensity 


All temporal fluctuations of variables for the reacting case are 


post-processed with limited number of sets of solutions stored 
due to lack of disk space. 
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Reacting case Nasal 


o—+ Exp. — KIES — K-LES-LEM — TENS 


i. 


* Experimental —— K-LES ~~ K-LES-LEM 


i 


0.35 0.70 0.35 0.7 0 0.35 0.7 0 0.35 0.7 0 0.35 0.7 
| a Ri” Aa 3 ae al x/D = 0.375 wD=0.95 xD=153 x/D=3.75 x/D =9.4 
wWD=1.53 x/D=3.75 


Transverse profiles of normalized mean Transverse profiles of normalized axial 
axial velocity at different axial locations. RMS velocity at different axial locations. 
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Reacting case 


—— K-LES-LEM — TENS 


* Experimental —— K-LES 


25 0 O25 OO 0.25 
x/D = 9.4 


04-04 0 0. 0 0.25 0 0.25 0 02 
WD=0.375 x/D=095 xD=153 xD=3.75 


Transverse profiles of normalized 
RMS Y-velocity at different axial 
locations 


04-04 0 04-04 0 04-04 0 


O04 0 
x/D = 3.75 w/D=9.4 


x/D = 0.375 x/D = 0.95 x/D = 1.53 
Transverse profiles of normalized mean Y- 
velocity at different axial locations 


SciTech 2017, 55'™ AIAA ASM Conference, Grapevine, TX. 01/12/2017 www.nasa.gOov 20 


National Aeronautics and Space Administration 


Reacting case 


— K-LES-LEM —— TENS 


y/D 
=) 


0.15 0 015-015 0 015-015 0 O15 -0.1 0 


0.1 -0. : 9 ) a 
WD = 0.375 x/D = 0.95 Say 184 wD =37 0 500 1000 1500 2000 0 500 1000 1500 20000 500 1000 1500 2000 


5 : x/D = 3.75 x/D = 8.75 x/D = 13.75 


3 


Transverse profiles of mean Reynolds stress Transverse profiles of mean temperature (K) 
at different axial locations from a two-step at different axial locations from a two-step 
global kinetic reaction model for propane . global kinetic reaction model for propane . 
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Concluding Remarks o 


¢ Anon-reacting case for the bluff-body flame holder configuration is 
investigated with TFNS approach on a 1600k polyhedral grid that 
refined from a 800k grid downloaded from the workshop webpage. 
The numerical solutions moderately predict the non-reacting flow 
field in terms of mean profiles. 

¢ For the reacting case, the bluff-body flame holder configuration is 
investigated with K-LES, K-LES-LEM and TFNS approaches on the 
original 800k grid downloaded from the workshop webpage. 

¢ All three reacting simulations predict that the flame is symmetric. All 
three reacting simulations produce time averaged axial velocities 
that are comparable to the experimental data in the recirculation 
zone. 

¢ The time averaged Y-velocities are predicted less accurately 
compared to the experimental data in the recirculation zone. 
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Concluding Remarks o 


¢ All three reacting simulations over-predict the mean temperature 
due to global mechanism. All reacting turbulence fluctuations are 
poorly predicted due to not sufficient data saved. 


¢ Finally, for the current work, the most important impact that 
determines the quality of the reacting simulations is the lack of 
the proper reaction kinetics. 

¢ For the next workshop, the coding effort to store unsteady data 
for many specified profiles inside the CFD code is pending. 
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